NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



NASA 

Technical Memorandum 80656 

(HlSA-TH-80656) SAIN BATE BANGE PROFILING N80-25989 

FROM A SPACEBOBNE RADAfi (NASA) 28 p 

HC AO 3/11 F AO 1 CSCL 04B 

Unclas 

G3/47 22944 

Rain Rate Range Profiling 
from a Spaceborne Radar 


Robert Meneghini 


MARCH 1980 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 



TM 80656 


RAIN RATE RANGE PROFILING FROM A SPACHBORNK RADAR 

Robert Mencghini 


MARCH 1^80 


GODDARD SPACE E'EIGIIT ( I NTER 


Greenbelt. Maryland 20771 
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Robert Meneghini 

ABSTRACT 

At certain frequencies and incidence angles the relative invariance of the surface scattering prop- 
t ties overland can be used to estimate the total attenuation and the integrated rain rate from a space- 
borne attenuating-wavelength radar. In this paper the technique is generalized so that rain rate pro- 
files along the radar beam can be estimated, i.e. rain rate determination at each range bin. This is 
done by modifying the standard algorithm for an attenuating-wavelength radar to include in it the 
measurement of the total attenuation. Simple error analyses of the estimates show that this type of 
profiling is possible if the total attenuation can be measured with a modest degree of accuracy. 


RAIN RATE RANGE PROFILING FROM A SPACEBORNE RADAR 

Robert Meneghini 

Introduction 

To achieve adequate resolution in remotely sensing precipitation from space, investigators have 
reconsidered the potential of attenuating wavelength radars. As first shown by H' ichfeld and Bordan 
[ 1 ] , the presence of attenuation can severely limit the accuracy in the determination of rain rate. To 
circumvent this problem, a number of dual measurement techniques have been proposed. Suitable 
forspacebome radars are the dual-attenuating wavelength [2,3) and the surface reference techniques 
[41 . For ground based measurements, additional information is provided by an independent meas- 
urement of rainfall rate at some point along the radar beam [ 1 ] or by an estimate of attenuation be- 
tween the radar and a fixed target [51 . One other method, the combined radar-radiometer technique 
[6J , appears to have applications to both ground based and spaceborne measurements. 

The fixed target and surface reference methods were originally analyzed from the perspective of 
providing average or integrated rain rate [4,5] . The observation upon which the techniques rest is 
that the attenuation factor can be found by forming the ratio of return powers from the target in the 
presence and absence of precipitation along the beam. The estimate of integrated rain rate then fol- 
lows upon using an empirical law that relates rain rate to attenuation. 

In these techniques, however, no use is made of the information provided by the measured re- 
flectivity factors. The major purpose of the paper is to show that these reflectivity factor measure- 
ments along with the total attenuation can be used to give range profiled rain rates, i.e. rain rate de- 
termination at each range bin. 

To do this, we can begin with either the Hitschfeld-Bordan [ 1 ] or the iterative algorithms of 
Ormsby [71, both of which depend only on the measured reflectivity factors and several empirically 
derived parameters. Although these algorithms themselves provide bin by bin estimates of rainfall 
rate, they are susceptible to large errors whenever the total attenuation is significant [ 1 , 7-101. For 
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spaceborne radars, the resolution requirement makes it desirable to operate at higher frequencies 
where the attenuation will be non-negligible for all but very small rain rates. Consequently these 
techniques in themselves are of limited applicability from space. We can, however, reformulate them 
so as to account for the fact that the total attenuation has been estimated by means of the surface 
reference method mentioned above. By incorporating this additional information into the algorithm, 
we can usually reduce the errors in the estimated rain rate. 

The foregoing correction procedure is reminiscent of that proposed by Hitschfeld and Bordan [ 1 ] . 
In their technique an independent point measurement of rain rate along the radar beam (obtained, for 
example, by means of a rain gage) is used to correct the rain rate estimates up to that point. The dif- 
ference between the two procedures is that one uses a point measurement of rain rate, the other an 
estimate of total attenuation It should be noted that Hitschfeld-Bordan correction technique, like 
the fixed reference target technique, applies only to ground based radars pointed near the horizon. 

The surface reference technique which uses the ground as a reference target is suitable for spaceborne 
or airborne radars. 

In this paper we will consider only those estimates which are based in the Hitschfeld-Bordan 
algorithm since the algorithms derived from it have been shown to be superior to those based on the 
iterative algorithms. Five estimates will be analyzed: the original Hitschfeld-Bordan (H-B) algorithm, 
the H-B algorithm that incorporates an estimate of the attenuation, and the H-B algorithm that in- 
corporates a point measurement of rain rate. For each of the last two cases, two distinct estimates 

are derived depending on whether the additional information (attenuation or point rain rate) is used 

% 

to correct, the calibration constant or, the parameter a where k = alP , k is the attenuation coef- 
ficient in db/km and Z is the reflectivity factor in mm 6 /m 3 . 

Once the five estimates have been written down a simple two-part error analysis is carried out. 

In the first part we assume that the errors in certain parameters can be characterized by the true value 
plus a constant offset term. This assumption is helpful in isolating the effects and assessing the rela- 
tive importance of the various errors. In the second section, the parameters are assumed to be random 
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with mean values that may be offset from the true value. With these less stringent assumptions, the 
error analysis should provide an indication of the conditions under which the estimates can be used. 

Rain Rate Algorithms 

For the rain rate in the j th range bin, the Hitschfeld-Bordan estimate can he written [ 10] , 




where 


j “1 ” b /0 
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K = 0.2 a fin (10) 
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The measured reflectivity factor, Z m ., the unattenuated reflectivity factor, Zj, the radar return power. 
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P R ,. and the calibration constant C are related via the equations, 


\ s A i h ‘ r i V c «> 

where Tj is the distance from the radar to the j th range bin, and where the attenuation factor at the 
j ,h range bin, A^, is defined by 

Aj = 10 k(xWx 


(3) 


where k is tlie attenuation coefficient in db/km. Using the empirical relation k = a7?, (3) can be 
written as 


A, = 10-°- 2a jT rj 


(4) 


Equation ( 1 ) follows from the Z-R relationship 


or, on using (2), 


Pj = aZj* 


*1 * »Z* /A? 


(5) 




and die fact that 


Af = 


j 

i- 1 


(6) 


which is a consequence of (2) and (4). The quantity s in ( 1 ) and (6) is the radar range resolution 
which enters through approximating the integral in (4) by a summation, i.e. 
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Assume now that a measurement of the attenuation factor out to the n ,h range bin has been ob- 
tained; that is. we have available an estimate A n for the quantity defined by (3). There are several 
ways to incorporate this additional information into (1 ). We discuss two possibilities. 

Case I. To account for an offset in the calibration constant C we multiply the set {Z,^} in (1) by a 
constant A. Notice that this is equivalent to multiplying C by A" 1 in (2). The Hitschfeld-Bordan 
estimate then becomes 

: -b/0 


R HB =a(AZ m .) b 
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Identifying AZ m . as the corrected form of Z m ., then from (5) and (6) 
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Letting j go to n in the above formula, using for A n the measured value A n , and solving for A gives 

I " | 1 10 

A = 


(1-AJ) 

i 

This quantity is inserted back into (7) which yields a new estimate of rain rate that we label with a 
superscript 1 ; 

r; .;< j {K ( is|„-A;r'S n -S 1 l}-''" j< n (8) 
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where 
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Equation (8) is an alternate form of the Hitschfeld-Bordan algorithm in which the measured quantity 
Ag has been used to correct for errors in the calibration constant. 

Case 2. To compensate for offsets in a, we replace K by KA in (1 ) and identify the corrected form 
of K by the latter quantity. Using (6) with K -* KA and proceeding in the same way as Case 1 , then 

A = (l-A')/K0sS n 

so that the second estimate, labeled with a superscript 2, becomes 

Rj < 2) * a z ^.j 1 1 - (1 - Aj|) Sj/S n ) j < n (10) 

As there is only one additional measurement (that of A,,) we can attempt to eliminate offsets in 
either a or C but not both. Within the approximations of the error analysis presented later in the 
paper, the rain rate as determined from (10) is independent of offsets in a but dependent on calibra- 
tion constant offsets. Exactly the opposite is true of (8). 

In their original paper [ 1 ] , Hitschfeld and Bordan proposed a correction scheme via an inde- 
pendent measurement of rain rate along the radar beam. This may be obtained, for example, by 
means of a rain gage located under the radar beam. We denote this measured rain rate at or near the 
n th bin by R ( . Again we consider two cases. Multiplying the set {Z m J in (1 ) by A, then 
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Letting j increase to n and using R f for R n then 

S,-.(AZ» / 
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Solving ( 1 2) for A and substituting this into ( 1 1 ) gives (Case 1 ) 

Rj 3) * aZ^ R f Z* mn +W Rj^ (S n - S,)) - b/ * j < n (13) 

Alternatively, K can be replaced by AK and an identical procedure carried out, yielding (Case 2) 
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(14) 


R| 4) = aZ^ [ 1 - 7 I ~ b/(} j < n 

with 

TMl-taZ^/R^HSj/SJ (15) 

Error Analysis 

Tlie sources of error in the above estimates arise from random fluctuations or offsets in the 
reference measurement of rain rate or attenuation, in the parameters appearing in the k-Z, Z-R rela- 
tions, and in the radar calibration constant. In the empirical k*Z, Z-R relations, the major source of 
error is caused by the dependence of scattering and attenuation on the temperatu-e and drop size 
distribution (DSD) of the rain. Atlas and Ulbrich |3] have noted that the DSD exhibits large syste- 
matic changes with storm type imposed upon a random component representing the spatial and 
temporal fluctuations within each class of precipitation. Thus, in general, the assumed parameters in 
the k-Z, Z-R relations will differ from the true values by an offset plus some random component of 
zero mean. 

Sources of error that affect the accuracy in determining Z m are the offset in the calibration 
constant, the sampling errors caused by the random nature of the scatterers, and the receiver noise. 

In the simple error analysis presented here, sampling errors and finite signal to noise ratios will be 
neglected; i.e. we assume that the radar has been designed and is operated under conditions for which 
these errors are small. 

Another type of error arises from the uncertainties in the reference measurement. In the fixed 
target and surface reference methods, the measurement of the attenuation factor will be corrupted 
by the presence of precipitation in the range bin containing the target. If a rain gage is used for an 
independent determination of rain rate errors can be caused by a lack of spatial averaging and the 
existence of too long a temporal average. Additional errors are incurred because of the difference in 
rain rates between the surface and the radar beam. 

The estimates given by (8) and ( 1 0) depend on the attenuation factor and the n reflectivity 
factors. As mentioned above, these quantities can be obtained if either the surface reference or fixed 
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target method is used. If, however, a radiometric measurement of brightness temperature were avail- 
able, the attenuation factor could be estimated from this. The radar then would be used only to 
supply the reflectivity factors. 

For spacebome application over ocean the advantage of this more complicated technique is that 
such a radar-radiometric sensor could provide rain rate profiling without recourse to the surface 
reference technique - a technique that is expected to be less accurate over ocean than over land. 

The errors in the radiometric measurement of precipitation, like those of the radar, depend on 
the differences between the actual and assumed values of the temperature and DSD of the scatterers. 
For spacebome applications the emission and scattering of radiation from the earth’s surface contribute 
to the observed brightness temperature with the consequence that quantitative measurements of rain 
rate can be made only over cold backgrounds, e.g. oceans. One disadvantage of the radiometer is its 
inability to measure the storm height, a quantity needed to obtain the attenuation factor. In a radar- 
radiometric sensor, the radar would be used to estimate this quantity. Discussions of these and other 
sources of radiometric error are given in [11-13]. 

As the attenuation factor may be obtained by several methods, each subject to different sources 
of error, it would seem necessary to analyze each in order to characterize the errors in A. Previous 
investigations [4,5] have shown that fairly accurate estimates of A can be made by means of the sur- 
face reference or the fixed target techniques. We therefore choose the offset and standard deviation 
of A to be typical of the values appropriate to a spacebome radar using the surface reference tech- 
nique. The results to be presented can also be viewed as determining the accuracy to which A (or R g ) 
must be known in order that the profiled rain rate be within the given error bounds. 

To analyze the estimates it is necessary to characterize the errors in the k-Z, Z-R relations as well 
as the measurement errors in A and R ( . It is convenient to distinguish two cases where the data set 
consists of measurements for which: the DSD is relatively constant or the DSD is variable. In the 
former case, the parameters in, say, the k-Z relation can be written in the following way (where 
k * aZ 9 ): fixing 0, we define a best fit value Oj. that can be determined from the true DSD. The dif- 
ference between a T and the assumed value a is a fixed offset error, E a 
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or 


a,. * a + E a 


a » 5 Q a T (16) 

with 

6 a = (a T -E a )/a T (17) 

When the data consists of measurements over many storms or over a storm in which the DSD changes 
in time, 6 0 should be interpreted as a random variable. For measurements taken over a short period 
of time (one record) we assume that there is no range variation in a T . The variability in Oj. then 
arises only from the fluctuations in the mean DSD from record to record. With this assumption we 
can still use (16) with 6 a taken to be random. 

Interpreting the quantities a, A n ami R g in a similar manner we can write 

a = 5 a a T (18) 

< 19 > 

R f * S R R, = 5 r R„ (20) 

where R n , A n are the true rain rate and the true attenuation factor at the n ,h range bin. As with 5 a , 
the quantities 6 a , 6 A . 6 R will he chosen either as constants or random variables depending on the 
nature of the data set. 

The radar calibration constant can be written in the form C = 6 f Cj where 6 C is a non-random 
quantity. Using this we can write at the j th range bin 

■ Z„,0 l'«r < 21 > 

where Z (T) denotes the true measured reflectivity factor. Notice that as the receiver noise and the 
sampling errors have been neglected. I is subject only to errors in the radar calibration constant. 

Mlj 

We now use the above equation to write the estimates as function of the five quantities (fi a , 6 a , 
5 a , 5 r , 6 f ). To simplify the results, the summations over Z [n <T) are eliminated by means of the 
relation 

j 

Af =1-0.2 («n 10) a T /.is ) Z 0 m . (T) 
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where Aj is the true attenuation factor up to the j* 11 rente bin. This equation can be found from (6), 
(17) and (19) and by recognizing that the Hitschfeld-Bordon estimate is exact in the absence of 
errors. We also use the fact that (5) 

*1 •*!*», (TV'S* 

where all the above quantities are the true values. 

Omitting the simple computations and using the same notation for the estimates as above the 
ratios of the estimated rain rates to the true rain rate at the j ,h range bin are given by 

5, 

Rj H ' b /Rj* — 11 + e ]‘ b/tf j < n (22) 

e 0 «Aj* <1-Af)(l -*./*£) 

R(l)/I S* ~ n-e,!- 1 j<n (23) 

O 

«i a«> 


R®/R, - "7“ 


Af(I-A») 


b If 


R| ,, /R j * a. « R 


C (6^ - 1 ) + A? (1 - 6® Aj)J 

b It 


j<n 


6 J /b C 3 + 6 a 6 R b 


j<n 


R n“*R 


(24) 


(25) 


e, « A*/(Af - A?) 


k l 4)/R rzr 


'i 


Af (1 -AS) 


“1 b If 


(Af-AS) + e 4 Aj(l-Af) 


j<n 


(26) 
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Certain features of the above estimates are readily deduced. For the Hitschfeld-Bordan estimate 
(22) we summarize some of the observations made by previous investigators. From the definition of 
the attenuation factor, it is apparent that the quantity A^. 0> 0. is an exponentially increasing func- 
tion of rain rate and range (i.c. of j). As Aj -0 increases, the absolute value of e 0 in (22) will also in- 
crease. If 6 a /6£ < 1 . the sign of e o will be positive and Rj* “ B will tend to be less than the true rain 
rate, R^. If 5 a /5^. > 1. e o will be negative and RM' b will generally be larger than Rj. Notice that as 
e o -*-l the estimate fails to exist. Ir. general, the presence of the A"** factor shows that whenever 
the total attenuation up to the j ,h range bin is significant, even modest errors in o or C can result in 
large errors in Rj* ' B . 

For the estimates Rj 1 \ Rj 2 * nearly the opposite situation exists: for small values of attenuation 
(out to the n ,h range bin), the estimates arc subject to large errors which then decrease as the attenu- 
ation becomes larger. For the R^ estimate (23) it is easily shown that: 

I < 1 1 -c, I ' 1 < (1 - 7)' 1 for 6 a < 1 (27) 

(1 + I7I)" 1 < 1 1 -e, I " 1 < 1 for « A > 1 (28) 

with 

■> - 1 1 - «/< I - JS- AS » 

We impose the condition 8^ A^ --- < I which is equivalent to stating that only positive values of 

measured attenuation are admissible. 

When the attenuation is large, >> 1 . then 7 ~ I - 5 A so that 1 1 - 6 , | _I is bounded between 
I and S”** for fi A < 1 and between 5*" and 1 for 6 A > 1 . Thus, for large values of attenuation, the 
maximum error in the term 1 1 - c, | is nearly proportional to the measurement error in the atten- 
uation factor. A n . If the attenuation is small, however. < 27) and (28) show that the error in [I -e,l 
and thus in Rj 1 * can be large. IT s is same general behavior is true of the Rj 2) estimate. 

From these considerations and anticipating the numerical results, we conclude that whenever 
the attenuation is significant the R (l) , R 1 -* estimates are generally superior to the R HB estimate. 
For very small values of attenuation, where R‘ ll . R *’ 1 are subject to large errors, the R H B esti- 
mate is preferable and often yields satisfactory results. 
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The estimates ffc (3) , ft (4) differ from the others in the sense that the rain rate at the n th bin has 
been controlled rather than the total attenuation. From (25) the following approximation can be 
deduced, 

fi./*!* e 3 « > 

Rj 3 ) /Rj « ■< 

[^r e 3 » 1 

The quantity e 3 is small when the total attenuation is large and j is not close to n;e 3 is large 
when the total attenuation is small or j is near n. For other values of e 3 , Rp J /Rj will be intermediate 
to these extremes. Two characteristics noteworthy of the R (3) and R (4) are that the errors are 
bounded both for large and small values of the attenuation and that for e 3 » 1 the estimate is 
independent of errors in the quantity ‘a’, a term which can make a sizeable contribution to the total 
error. 

Results 

Figures 1 through 4 illustrate the behavior of the estimates as a function of range for selected 
values of rain rate and offsets. The radar range resolution is chosen to be equal to .25 km with n = 20 
so that the distance from the radar to the n th range bin is 5 km. In all the figures the five ratios Rj/Rj 
of (22) - (26) have been plotted as functions of radar range. Notice that since Rj is the true rain race, 
an ^..crless estimate would give a ratio of 1. The wavelength, chosen equal to 0.86 cm, enters pri- 
marily in the determination of the attenuation: the results obtained for a 5mm/hr rain rate at X = 0.86 
cm will be nearly the same as those for R = 8 mm/hr at X = 1 .24 or R = 65 mm/hr at X = 3.2 cm since 
the attenuation is comparable in all three cases. 

In figure 1 , the rain rate is fixed at 5 mm/hr with e A = e a = e R = -0.25, e t = 0, e c = .25 where 
e denotes the relative error in the quantity appearing in the subscript, e. g. e A = 1 - 5 A = (A n - A n )/A n . 

The ratio 5 a /8£ with ]3 = .905 is significantly greater than 1 and consequently the R H ' B becomes 
much greater than the true rain rate within a short distance from the radar. The estimates R (1) through 

R (4) exhibit more gradual changes in range and much smaller errors than in R H ' B . 
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Figure 3. R*/R versus range: (R 1 = i ,h rain rate estimate, R = true rain rate) 
R = 5 mm/hr; X= 0.86 cm,(5 i( 5 A , & a , 6 C ,5 R ) = (1, .75, .75, .75, .75) 


RANGE (KM)- 


Figure 4. R*/R versus range; (R 1 = i ,h rain rate estimate, R = true rain rate) 
R = 20 mm/hr; X = 0.86 cm, (5^ 5 A , 5 a , 5 C , 6 R )«(|, .75, .75, .75, .75) 


In figures 2-4. the relative errors are chosen to be e A = e a = e R = e c = 0.25, e a = 0 so that 
8 0 /8£ ~ -97. In figure 2. assuming R = 1 mm/hr, the errors in R (2) and especially in R (1) are larger 
than those of R H ' B , R l3) or R (4) . 


In figure (3), with R = 5 mm/hr, R**‘ B decreases more rapidly with range but as 5 a /6£, is close 
to 1 , the overall error remains small. As the rain rate is once more increased to 20 mm/hr (figuie 4) 
the errors in R l 1 * and R (2 > are considerably reduced whereas R H " B is rendered useless because of high 
attenuation. 


The dependence ot the estimates on rain rate, as shown in figures 2-4, is general: for low values 
of attenuation lv*' B is generally less in error than R (1 * and R* 2 * while the opposite is true at moder- 
ate to high values of attenuation. Hie transition point, however, is a sensitive function of the ratio 
funeral as the derivation ot 5 (i / 5j*. from 1 increases, IT 1 ' and R (2 * will become preferable 
to R» « at progressively smaller values of rain rate. In the foregoing discussion, the values in a, a, 

A n were assumed to be equal to the best tit value plus a constant offset. As stated previously, 
this is a reasonable assumption when the DSD is relatively constant over the data set. Of greater 
interest is the situation in which the DSD varies randomly between the data records. In this case an 
interpretation of 8 (i . 6 A , 5 R as random variables is more appropriate. Under this assumption the 
Rj/Rj ratios are themselves random. At X - 0.86 cm we can assume without much error that b = p = 1. 
Under this assumption, the mean and standard deviation of each of the rain rate estimates can be 
found by means of straightforward integrations. For simplicity we will present the statistics only for 
the R ,,B . R (1) , R 1 ' estimates. The analytic expression for the mean and variance of these quantities 
are given in the Appendix. 


Once the mean and standard deviation of A , a and « are chosen, we compute the absolute value 


ot the mean and the standard deviation of the Rj/R ratios as a function of range for a fixed, uniform 
rain rate, i.e. Rj ; R. To simphly the presentation of the results, an additional average is performed 
over range bins j - 1 to n: that is, we compute the quantities \f ( { K ), (7(R) where 
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R' 


i i \ 
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where the angular brackets denote the expectation. The notation Rj represents the i th rain rate esti- 
mate at the j lh range bin where R 1 denotes either the R H * B , the R (1 * or the R (2 * estimate. The cri- 
teria for a good estimate is that 5 l be small and that be close to unity. 

We again choose X = 0.86 cm, r, 4 = 5 km and s = 0.25 km so that n = 20. In tables 1 and 2, 

(Mj . 5j), (M 2 , ff 2 ) and (M HB , 5 HB ) are listed as functions of rain rate for the three cases 6 C = (1, 
1.25, .75). These values correspond to calibration errors of 0 db, .97 db and -1.25 db respectively. 

In table 1, o g = o a = o A =0.125 where o g is the standard deviation of o g , etc. In table 2, o g = o a = 
o A = 0.25. For both sets of calculations the mean values of 6 g , 6 a and 5 A are assumed to be un- 
biased. Writing <6 g > = 1 - b g , etc. The above condition is equivalent to setting b g = b a = b A = 0. 

Since the Rj 1 * estimate is independent of calibration errors, M,, are independent of 6 C . 
Therefore, in tables 1 and 2 it suffices to list Mj, 5j for the 6 C = 1 case alone. For large rain rates, 
the mean or standard deviation of the Hitschfeld-Bordon estimate might not exist at one or more 
range bins. In these instances the range averaged statistics are undefined. We represent this in the 
tables by writing dashes for the M H . B , o HB . 

The results show that ft H ' B is useful only at small values of attenuation, i.e. low rain rates. Of 
the three estimates llj 2 * has, in general, the smallest variability about its mean value. Relative to the 
Rj 1 * estimate, its disadvantage is the dependence on the radar calibration error. For example, at 
8 C = 1.25, the mean value of R (2 * is approximately four-fifths of the true rain rate while at 5 C = 0.75 
it is approximately one-third greater than the true rain rate. In tables 1 and 2 the R (1 * estimate is 
nearly unbiased for the rain rates considered. Notice that the increase in o g , o a , o A from table 1 to 
table 2 has only a slight effect on the mean values but results in a significant increase in the variance 
of all three estimates. 

In table 3 we have set o g = o a = o A = 0.25, 6 C = 1 and b g = 0. As before (Mj, Oj), (M 2 , 5 2 ) an d 
(M h . b , a H . B ) are given as functions of the true rain rate. In the first example, b a = b A = 0.2 so that 
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Tabic 1 ^ ^ ^ 

Range Averaged Mean and Standard Deviation of R (1) , R (2) , ^ H ' B 
versus the True Rain Rate 
°«* °a> °a = 0-125 
b ,,b a ,b A =0 


Rain Rate 

5 C = 1 

(mm/hr) 

MpSj 

M 2 ,o 2 

MH-B* ff H-B 

1 

1.03, .388 

1.0, .144 

1.0, .133 

2 

1.02, .252 

1.0, .142 

1.01, .168 

3 

1.02, .238 

1.0, .14 

1.06, .269 

4 

1.02, .211 

1.0, .138 

— — 

5 

1.02, .206 

1.0, .137 

— — 

10 

1.02, .197 

1.0, .133 

• - - 

15 

1.01, .25 

1.0, .141 

— — 

20 

.976, .339 

.99, .179 

- - 


«c = 

1.25 


1 


.804, .115 

.756, .098 

2 


.803, .113 

.707, .101 

3 


.803, .112 

.654, .106 

4 


.803, .111 

.599, .112 

5 


.802, .110 

.545, .115 

10 


.802, .106 

.336, .095 

15 


.8, .113 

.229, .067 

20 


.79, .143 

.172, .05 


5 C = 

0.75 


1 


1.34, .192 

1.5, .217 

2 


1.34, .189 

— - 

3 


1.34, .187 

— - 

4 


1.34, .184 

— - 

5 


1.34, .183 

— - 

10 


1.34, .177 

- 

15 


1.33, .188 

— 

20 


1.32, .238 

- 
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Table 2 ^ ^ ^ 

Range Averaged Mean and Standard Deviation of R (1 *, R (2) , R H ‘ B 
versus the True Rain Rate 


<VO«,o a -°.25 

b .,b«,b A -0 


Rain Rate 
(mm/hr) 

M l.*l 

M 2 ,5 2 

m h-b* 5 h-b 


«C 

® 1 


1 

1.15, .78 

1.02, .297 

1.01, .269 

2 

1.11, .544 

1.02, .292 

1.06, .38 

3 

1.1, .482 

1.02, .287 

— — 

4 

1.1, .456 

1.01, .283 

— — 

5 

1.09, .446 

1.01, .279 

— — 

10 

1.08, .427 

1.01, .269 

— — 

15 

1.08, .436 

1.01, .266 

— — 

20 

1.04, .538 

.996, .286 

_ _ 


6 c = 

1.25 


1 


.817, .237 

.759, .196 

2 


.815, .233 

.718, .207 

3 


.813, .229 

.682, 236 

4 


.812, .226 

.667, .33 

5 


.811, .223 

— — 

10 


.807, .215 

— 

15 


.805, .213 


20 


.797, .228 



«c = 

0.75 


1 


1.36, .395 

1.52, .449 

2 


1.36, .389 

- 

3 


1.36, .382 

— — 

4 


1.35, .377 

— — 

5 


1.35, .372 

— -- 

10 


1.34, .359 

— 

15 


1.34, .355 

— _ 

20 


1.33, .381 

- - 
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Table3 

Range Averaged Mean and Standard Deviation of R (1 \ R (2) , R H ' B 
versus the True Rain Rate 


V °a * 0 25 
*C " ! > b . "° 


Rain Rate 
(mm/hr) 

Mj, 8 

l 

M 2 , 

5 2 

m h-b* 5 h-b 



b a = b 

>A - 0-2 



1 

2.2, 

1.5 

1.15, 

.356 

.951, .250 

2 

1.81, 

1.0 

1.13, 

.345 

.908, .282 

3 

1.67, 

.865 

1.12, 

.334 

.891, .379 

4 

1.61, 

.791 

1.11, 

.325 

- - 

5 

1.57, 

.765 

1.10, 

.318 

- - 

10 

1.49, 

.706 

1.06, 

.295 

- - 

15 

1.46, 

.703 

1.05, 

.286 

- - 

20 

1.40, 

.787 

1.03, 

.303 

- - 



b a = t 

, A =-0.2 



1 

.589, 

.484 

.93, 

.26 

1.08, .292 

2 

.730, 

.331 

.93, 

.258 

- - 

3 

.777, 

.304 

•94, 

.257 

- - 

4 

.80, 

.296 

.94, 

.256 

- - 

5 

.814, 

.297 

.95, 

.255 

- - 

10 

.842, 

.298 

.966, 

.253 

- - 

15 

.85, 

.311 

.975, 

.254 

- - 

20 

.825, 

.397 

.972, 

.276 

- - 
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<8 # > * <$ A > * 0.8; in the second example b 0 ■ b A ■ -0.2 so that <8 a > * <8 A > ■ 1 .2. In both cases, 
the R (1) estimate is significantly offset from the true rain rate; moreover, in the case of b A ■ b 0 * 0.2 
the standard deviation is very large when the rain rate is small. 

For both examples, the R (2) estimate is clearly superior to the other two. We conclude that when 
the radar is well calibrated, R (2) yields the most accurate results of the three estimates considered. 

Summary 

Several recent papers have shown that integrated rain ra ts estimates can be obtained by means 
of an attenuating wavelength radar. These estimates arc found by means of a measurement of the 
attenuation factor, A n . In this paper we have shown that range profiled rain rate estimates also can 
be deduced by using A,, along with the Hitschfeld-Bordon algorithm and the measured reflectivity 
factors {Z m j}- In the derivation it was seen that two estimates, R (1) , R (2) could be obtained from the 
measured data, the first of which was independent of erron in the calibration constant, C, the latter 
of which was independent of errors in a, where k - aZ 9 . If an independent measurement of rain rate, 
Rj, is made at some point along the radar beam two additional estimates of rain rate, R (3) , R (4) , 
analogous to R (1 \ R (2) , can be deduced. 

This latter method is essentially the correction scheme proposed by Hitschfeld and Bordan. 

For purposes of comparison we added to these four estimates the original range profiled estimate 
derived by Hitschfeld and Bordan, i.e. R H ’ B . Since R H ' B is a function only of {Z ml } this estimate 
is the only one of the five that can be employed if neither A n nor R g are available. 

To investigate the qualitative behavior of the estimates, a simple error analysis was performed. 

The major findings were: 

1 . When the total attenuation is signi^ant, the estimates that incorporate the additional 
measurement of either A„ or R,j are superior to the $ H * B estimate if these quantities 
can be measured with a modest degree of accuracy. 

2. The R (, \ R (4) estimates yield fairly accurate measurements of rain rate for both small and 
large values of total attenuation. On the other hand, the R u \ R <2) estimates often suffer 
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from large errors when the total attenuation is small. The drawback of the R (3) , R (4) esti- 
mates is that they can be applied only to ground based radars where an independent meas- 
ure of point rain rate can be obtained. 

3. The R (2) estimate is generally more accurate than the R (,) estimate when the errors in A n , 

C and a are of comparable magnitude. When the errors in C are much greater than those 
in a, the R-' • estimate is superior. 

Acknowledgments 

I wish to thank Dr. Jerome Eckerman for his assistance and encouragement. 

References 

[ 1 ] Hitschfeld, W. and J. Bordan, Errors inherent in the radar measurement of rainfall at attenu- 
ating wavelengths, J. Meteor., 11, 58-67, 1954. 

[2] Goldhirsh, J. and I. Katz, Estimation of rain drop size distribution using multiple wavelength 
radar systems, 9, 439-446, 1974. 

[31 Atlas, D. and C. W. Ulbrich, The physical basis for attenuation-rainfall relationships and the 
measurement of rainfall parameters by combined attenuation and radar methods, J. Rech. 
Atmos., 8, 275-298, 1974. 

[4J Eckerman, J., R. Meneghini and D. Atlas. Average rainfall determination from a scanning beam 
spaceborne meteorological radar, NASA Tech. Memo. 79664, Nov. 1978. 

[51 Atlas, D. and C. W. Ulbrich, Path and area-integrated rainfall measurements by microwave 
attenuation in the 1-3 cm band, J. Appl. Meteor., 16, 1322, 1331, 1977. 

[6) Hodge. D. B. and G. L. Austin. The comparison between radar- and radiometer-derived rain 
attenuation for earth-space links. Radio Science, 12, 733 740, 1977. 

[71 Ormsby, J. F. A., Some results on attenuation compensation analysis. Preprints 1 5 th Radar 
Met. Conf. Champaign-Urbana, 144-148, 1972. 

[8) Hildebrand, P. H., Iterative correction for attenuation of 5 cm radar in rain, J. Appl. Meteor., 
17, 508-514, 1978. 


20 


[9] Goldhirsh, J., Analysis of single and multiple wavelength radar systems for estimating slant 
path attenuation through rain and cloud at frequencies greater than 10GHz, APL/JHU, 
SIR77U-036, Oct. 1977. 

[10] Meneghini, R., Rain Rate estimates for an attenuating radar, Radio Science, 13, 459470, 
1978. 

[11] Wilheit, T. T., A. T. C. Chang, M. S. V. Rao, E. B. Rodgers and J. S. Theon, A satellite tech- 
nique for quantitatively mapping rainfall rates over the ocean, J. Appl. Meteor., 16, 551*560, 

1977. 

[ 1 2] Austin, P. M. and S. G. Geotis, Evaluation of the quality of precipitation data from a satellite- 
borne radiometer, Final Report, NASA-NSG-5024, Department of Meteorology, MIT, June 

1978. 

[13] Weinman, J. A. and P. J. Guetter, Determination of rainfall distribution from microwave 
radiation measured by the Nimbus-6 ESMR. J. Appl. Meteor., 16, 437442, 1977. 


21 


Appendix 

Statistics of the Rain Rate Estimates 

To calculate the mean and standard deviation of the normalized rain rate estimates Rj H /Rj, 
and il|' '/Rj. we must assume some knowledge of the random variables a, a, and A,,. As we 
are primarily interested in the qualitative behavior of the estimates, we can simplify the computation* 
by making the following assumptions. In die k-Z, Z-R relations, k ■ aZ f , R ■ aZ b , we take 1. 
Furthermore, we assume that the true values a T , Oj. and A n are uniformly distributed random vari- 
ables so that the probability density function of, for example, fi a can be written 


1 


2 2 -;(l-b l )-A # <« i <(l-b,) + A i 


W,)“ 


(Al) 


0; otherwise 
with similar expressions for 5 a , 5 A . 

We can rewrite (Al) by using the fact that A # * >/? o t , where o g is the standard deviation of 
6 g . The quantity b # in (A' ) o the relative bias in a, i.e. 


•••(V) 


where a T is the trw value, a is the assumed value and the angular brackets denote the expectation 
operator. Using A, 11 vT o > . ( Al ) can be written as 

1 


*«.)’ 


Id “b a ) VTo t ) < 5, < [(1 -bJ+v'J oj 


(A2) 


0; otherwise 


with 


<5 i )« 1 ~ b g * 1 -< 


f a T - a\ 


(A3) 


The mean and standard deviation c <f fue estimates can now be computed by means of iimple integra- 
tions. For the mean and mean square, of the Hitschfeld-Bordan estimate, we obtain 
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.* * * 


(A4) 


/gl\ ■ Ji^L h p-I&l 

\ «, / a^a^o-A,) U‘PoJ 

v/To, A,“» (1-Aj) 

P ° Sc+Aj^l-AjMVO-ba)) 

oj+d-b,) 2 

>1-3 

7j - 1 ♦ Aj* 1 (6 C -l) + Aj-» (1 - Aj) b # 
7 2 - A,- 1 (1 - A,) 



(A5) 


(A6) 

(A7) 

vA8) 


For the mean to exist, the condition -1 < p Q < 1 must be satisfied; for the mean square to exist 

». < -rf/3-rt 

For the estimate Rj'VRj the mean and mean square are 


with 



(l-b § ) 

UT7 


fin 


1 +P, 


1 -P, 


C 0+ C, fin 


l-q, 


1 +q, 


C 0 “ ( 1 — AJ" 1 )”* 

P, -N/Ta./d-bJ 

q, -\T5o a A n (1 - Aj‘* )/((! -A„) + b A A„ ( I - A’ 1 )| 
C, - Aj- 1 (I - A b )/2 v'3 o a A n (1 - A,' 1 ) 2 

0 3 2 +(»-b 1 ) 2 



(i b a ) 2 -3 a\ 


ic 2 *C, +C. fin 


> 

> 

to 2 3 

1 + q j 

> 


(A9) 


(A 10) 
(All) 
(AI2) 

(AI3) 


Id - A n )/A j (l -A' 1 )] 2 

M - 

2 1(1 -A B ) + d -A j ’ , )b A A n 1 2 - 3 ( o A (1 -Aj-bAJ 2 


(A14) 
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(AI5) 


Cj ■ (1 - \ n )/y/J o A Aj A„ (1 - Aj" 1 )* 

For both the mean and mean square to exist the following conditions must hold: 

•i<q, < i 


-l < p, < i 

a a <d-b a )/VT 


Finally, for the mean and mean square of Rj J, /Rj we have 

A|(l-A,)(l-b,) 


where 


fin 


*,/ ~ 2>/TA.(l-A,)« A « e 

Aj ( 1 ” A„)-b A \ a (I -Aj) 


L + i> 

i-p. 


where 


^J <2) \ 2 \ A / (1 " A n )2 (0 l + (I * b . )2) 

K R il) 

7 3 * Aj (1 " A n )- b A A n (1 - Aj) 


(AI6) 

(A 17) 
(AI8) 

CA19) 

(A20) 


For the mean and mean square to exist 

-1 < P 2 < 1 

°A < y 4 

The standard deviation of the above estimates is given by the square root of the difference between 
the mean square and square of the mean. 


